Enhanced plasticity in Zr-Cu-Ag-Al-Be bulk metallic glasses 
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Abstract 

The dependences of bending plasticity on compositional adjustment and structural 
relaxation in Zr-Cu-Ag-Al-Be bulk metallic glasses (BMGs) were systematically 
studied. We reveal that with the increase in the Zr content, the plasticity is enhanced 
at the expense of glass forming ability and the shear transformation zone (STZ) 
volume becomes larger, while structural relaxation results in the embrittlement 
accompanied with the further enlarged STZ volume. Average positron lifetime, 
representing average free volume, drops with the increase in the Zr content and 
structural relaxation. The results indicate that neither STZ volume nor average 
positron lifetime can solely describe the plastic deformability of BMGs. Analyzing 
positron annihilation lifetime spectra, a bimodal distribution of open volume size is 
observed. The smaller open volume shrinks with increased concentration and the 
larger open volume expands with decreased concentration due to the increase in the 
Zr content, while structural relaxation only reduces the size of smaller open volume. 
The effects of open volume size distribution on the STZ activation and shear band 
propagation have been discussed. It has been demonstrated that combining the STZ 
volume with open volume size distribution can explain the plasticity evolution with 
the composition and structural relaxation, which sheds light on understanding the 


deformation mechanism of BMGs. 


1. Introduction 


As a class of advanced materials, bulk metallic glasses (BMGs) are promising for 
future structural, chemical and magnetic applications due to high strength, excellent 
corrosion-resistance and good soft-magnetic property, as compared with their 
crystalline counterparts [1-3]. A number of BMGs with the critical diameters above 
20 mm have been reported in various alloy systems [4]. However, strain softening 
often causes catastrophic failure through unhindered shear banding, significantly 
limiting the actual application in structural materials [5-7]. Therefore, how to 
overcome the brittleness has become one of the most pressing issues in the 
community, and great efforts have recently been devoted to enhance the plasticity of 
BMGs. For instance, heterogeneous microstructures have been designed in virtue of 
the interaction between glassy matrix and second phase, e.g., crystalline particles [8] 
or amorphous phase [9], to restrict the propagation of highly localized primary shear 
bands, and tensile ductility of the composites can be achieved. By cold rolling [10,11] 
and surface shot-peening [12], the plasticity of deformed BMGs was also improved in 
subsequent mechanical tests, which was ascribed to the introduction of pre-existing 
shear bands and compressive residual stresses. 

In the monolithic BMGs, the plasticity normally correlates with the value of 
Poisson’s ratio v, i.e. larger v, allowing shear slip to occur readily before the structure 
collapse caused by normal stresses, results in better plasticity [13]. Lewandowski e al. 
[14] analyzed the plasticity and elastic constants of various BMGs, and found a 
critical value of v (0.31~0.32), at which the brittle-to-ductile transition occurs. With 
the increase in v above this critical value, BMGs become more ductile with the 
enhanced plasticity. In addition to v, the plasticity is also correlated with free volume. 
Chen et al. [15] found that faster cooling and minor Ti adding can enhance the 
compressive plasticity of Cu-Zr-Al BMG without obvious change in v, and attributed 
the improved deformability to higher free volume induced by high cooling rate and 
minor alloying. Such a free-volume effect on plasticity was proposed by Jiang et al. 
[16] to be caused by (i) introducing more sites favorable for the initiation or 
nucleation of new shear bands; (ii) easing the stress concentration in the areas by 


enhancing the atomic mobility and rearrangement. 


Tuning alloy composition has been demonstrated to be an effective way to 
significantly meliorate the plasticity of monolithic BMGs through adjusting v and the 
free-volume content. However, what are the intrinsic factors controlling the plasticity 
of BMGs is still unclear. When the applied stress exceeds the yield strength, the local 
rearrangement of atoms, which can accommodate shear plastic strain, takes place and 
eventually results in the formation of shear bands. Such a cooperative rearrangement 
of atomic clusters overcoming the saddle point of an energy barrier in potential 
energy landscape was termed as shear transformation zone (STZ) [5,17]. Recently, 
Johnson ef al. [18] proposed a cooperative shearing model (CSM) to interpret the 
plastic deformation of MGs, and the obtained STZ size is of importance for 
understanding the microscopic deformation mechanisms and macroscopic mechanical 
property. But the correlation between the STZ size and the plasticity in BMGs is still 
under debate. For example, Pan ef al. [19] measured the rate dependence of hardness 
through rate-jump nanoindentation testing to calculate the strain rate sensitivity, and 
characterized the STZ size based on the CSM theory. By comparing the STZ sizes of 
various BMGs, a good correlation between the STZ size and v exists, i.e. bigger STZ 
size corresponds to high v (better plasticity), and vice versa. However, Liu ef al. [20] 
measured the STZ sizes of BMGs by using dynamic mechanical analysis, and found 
that the value of v drops gradually with the increase in STZ size. Therefore, the 
controversial results on the dependence of STZ size on v in BMGs suggest that the 
details of the correlation between the plasticity and STZ size need to be further 
elucidated. 

Recently, the composition of ZraCu3014Ags36AlsBe7s (at.%), which can be cast 
into glassy rod with 73 mm in diameter by conventional copper-mould casting, was 
successfully developed, but this BMG still suffers from low plasticity [4]. The 
excellent glass forming ability (GFA) of Zra&Cuso.14A gs 36AlgBe7,s provides an enough 
wide window for the compositional optimization with good plasticity and high GFA 
simultaneously. In the present work, systematical composition tuning in the Zr-Cu- 
Ag-Al-Be alloy system around Zr46Cu30,14Ag8.36Al8Be7 5 was performed, and thermal, 
mechanical and density parameters of these newly developed Zr-Cu-Ag-Al-Be BMGs 
were presented. In order to further understand the microscopic mechanisms of various 
plastic behaviors, the STZ sizes of BMGs were measured by studying either the rate 


dependence of hardness during nanoindentation creep or the statistical distribution of 
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hardness. The correlation between the STZ size and the plasticity of BMGs with 
different compositions or structural states was carefully discussed. Meanwhile, the 


key role of free volume on the plasticity of BMGs was also pointed out. 


2. Experimental procedures 


The compositions of the alloys were chosen in the following way. As the pseudo- 
ternary composition map shown in Fig. 1, ZraCu3o14A gs 36AlsBe7: (hereafter referred 
to as Z0) with excellent GFA were selected as the starting point, and multiple lines 
according to certain chemical formula were drawn. The formula can be expressed as 
[Z0](oo-o(ZTiooyM,), (M=Cu, Al), y determines the position of the end point of the 
line, x indicates the different points along the line. The pre-alloyed ingots were 
prepared by arc-melting high purity (99.9%) Zr, Cu, Ag and Al with commercial 
Cuss 3Be44 7 alloy several times in a Ti-gettered purified argon atmosphere. Plates with 
the thickness of 1.5 mm and the width of 10 mm and rods with the diameter of 3 mm 
were fabricated by suction-casting into a water-cooled copper mold. Transverse cross 
section of the rod and top surface of the plate were examined on a Rigaku X-ray 
diffractometer (XRD) with Cu Ka radiation. Structural analyses of newly-developed 
Zr-Cu-Ag-Al-Be BMGs were also conducted by using synchrotron radiation x-ray 
diffraction (SR-XRD) at Shanghai Synchrotron Radiation Facility (SSRF). The beam 
size was 0.4x0.4 mm? and the wavelength used was 0.178324 A. High-resolution 
XRD patterns were recorded on an MAR 345 image plate with 8001200 um? pixel 
size. Scattering intensity /(q) was extracted by using the software package FIT2D 
[21]. Then, structure factor S(q), reduced radial distribution function G(r) and pair 
correlation function g(r) were obtained by PDFgetX2 according to the following 
equations [22]: 


G(r) = 4nr[(p(r) — po)] = SU 4-1) - sin(qr)dq (1) 
EP gI- G(r) 

g(r) =F si (2) 

S(q) = 135 f*r[g(r) - 1]: sin(qr)dr (3) 


where p(r) is the radial density function, pọ the average atomic number density, 
q = 4r sin(6)/A. Thermal analyses of the plate specimens were performed in 
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differential scanning calorimeters (DSC, Netzsch 404C and Perkin-Elmer Pyris 
Diamond) under a flow of purified argon at a heating rate of 20 K min”. 

Three-point bending specimens with dimension 20x4x0.8 mm” were cut from the 
as-cast plates, and polished by 1000-grit SiC paper and diamond polishing solution. 
Bending tests were performed on a universal testing machine (SANS CMT5205) at a 
displacement rate of 1.2 mm min” with the span of 10 mm. During testing, the load- 
deflection curves were recorded. For each composition, four bending tests were 
adopted to evaluate the plastic deformability. The tensile surface of the bent plate 
after fracture was observed by scanning electron microscope (SEM, HITACHI TM- 
1000). 

In order to structurally relax the as-cast specimen and reduce the free-volume 
content, annealing treatments at temperatures below the glass transition temperature 
(7¿) for 1 hour were conducted. After annealing, the specimens were water-quenched 
to remain in the glassy state. Positron annihilation lifetime spectrum (PALS) 
measurements of the as-cast and relaxed specimens were performed using a fast-slow 
coincidence ORTEC system with a time resolution of 190 ps (full width at half- 
maximum) in Institute of High Energy Physics, Chinese Academy of Sciences. Two 
identical pieces of 9x9x1 mm” were placed on either side of the Na positron source, 
and then this sample-source-sample sandwich was placed between the two detectors 
to acquire the lifetime spectra. The positron annihilation spectra were de-convoluted 
using the LT-9 software to set apart the physical meaning information, i.e., positron 
annihilation lifetime and intensity. The source correlations were included in the 
lifetime analysis. 

The acoustic longitudinal velocity (V7) and shear velocity (V;) of the cylindrical 
specimens (3 mm in diameter and 6 mm in height) were measured by an OLYMPUS 
5072PR model ultrasonic system with a measuring sensitivity of 0.5 ns. The carrying 
frequency of longitudinal wave and shear wave was 15 MHz and 10 MHz, 
respectively. The Young's modulus E, bulk modulus X, shear modulus G, and 
Poisson’s ratio v were derived from the acoustic velocities and the density data. Based 
on Archimedean principle, density measurements of the BMGs and corresponding 
crystalline counterparts were performed with a Mettler Toledo XS105 microbalance 
having a sensitivity of 0.01 mg. For the determination of the crystalline density, the 


as-cast specimens of different compositions with the diameter of 3 mm and the length 
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of 20 mm were sealed in an evacuated quartz capsule and annealed at 1000 K for 1 
hour, and then cooled in the furnace to make sure that the as-treated specimens were 
completely crystallized. 

Nanoindentation creep tests were carried out at room temperature on an Agilent 
Nano Indenter G200 fitted with a spherical indenter tip with the radius of 11.1 um, 
which was calibrated by measuring the elastic modulus of fused silica. Tests were 
conducted in load-control mode with the loading and unloading rates of 5 mN s”, and 
maximum load of 200 mN was adopted. The holding time at maximum load was set 
to be 10 s. For each composition, more than 400 indentations spaced 50 um apart on 
the surface were conducted. In nanoindentation experiments, thermal drift rate was set 
to be 0.05 nm s' to keep the reliability of the results, and the thermal drift correction 


was performed at about 10% of maximum load during unloading. 
3. Results and discussion 


3.1 Composition tuning towards good plasticity and its correlation with Poisson’s 


ratio 

Along the lines starting from Z0 in the pseudo-ternary composition map shown in 
Fig. 1, dozens of compositions were chosen to examine their GFA and mechanical 
properties. Fig. 2a shows the XRD patterns of the as-cast (Z0)(100-x)(Z1 100M) (M=Cu, 
Al) plate specimens with the thickness of 1.5 mm. With the increase in the fraction (x) 
of Zrioo,,M,, the GFA gradually decreases, manifested in the detected diffraction 
peaks from crystalline phases. Since the crystalline/amorphous composites was not 
the focus of this study, only seven compositions without crystals in the plate 
specimens were tested by DSC, as shown in Fig. 2b. For the sake of evaluating the 
GFA of these BMGs, the parameters of 7,5 (=7,/7;) [23] and y [24] defined as 
T,/(Tg+T)), which are regarded as indicators of GFA, were calculated and listed in 
Table 1 by using the characteristic temperatures measured by DSC, where 7; is the 
onset temperature of crystallization and 7; is the liquidus temperature. It is clear that 
the values of T, and y gradually reduce with the increase in the fraction of Zrioo.yM,, 
indicating the decreased GFA. For instance, the critical diameter for glass formation 
in (Z0)zo(Zr;sCuas)3o was found to be 7 mm (XRD result not shown here), much lower 


than 73 mm for Z0. 


Mechanical properties of the newly developed Zr-Cu-Ag-Al-Be BMGs were 
investigated through three-point bending tests. The plastic deflection value was 
adopted to characterize the deformability. Fig. 3a shows the flexural strength as a 
function of deflection for the as-cast specimens with different compositions. It is 
evident that Z0 does not exhibit any plastic strain, failing in a fully brittle way during 
bending. With increasing the fraction of ZriooyM,, the plasticity is remarkably 
improved, especially for (Z0)es(Zrs;Cus3)ss and (Z0);o(Zr;zsCuos)so which have the 
largest plastic deflection value about 1.2~1.3 mm. Fig. 3b-e depicts the SEM images 
of the as-cast Z0, (Z0)ss(Zr)is, (Z0)sz(Zrs;Cuas)ss and (Z0)79(Z1r7sCurs)30 specimens 


after fracture in bending, respectively. Only a few shear bands appeared on the tensile 
surface for the Z0 specimen, while slightly more shear bands were observed in 
(ZO)ss(Zr)ıs. Much higher density of shear bands can be clearly detected in the 
(Z0)es(Zre;Cus3)ss and (Z0)70(Zr7sCuzs)30 specimens, as compared with those in Z0 
and (Z0)ss(Zr)is, indicating the largely enhanced plasticity. Recently, Schroers et al. 
[13] developed a Pts; :Cui47Nis3P2»5 BMG with pronounced global plasticity under 
both bending and unconfined compression conditions, and attributed the formation of 
multiple shear bands to a high v of 0.42, which causes the tip of a shear band to 
extend rather than initiate a crack. The possibility of high v resulting in the improved 
plasticity was also investigated in the newly developed Zr-Cu-Ag-Al-Be alloy system. 
Ultrasonic measurements were carried out to obtain acoustic velocities for the as-cast 
specimens, and the calculated elastic constants were listed in Table 1. With the 
increase in the fraction of Zrioo,, M,, the value of v gradually augments. Among the 
compositions investigated, (Z0)ss(Zrs;Cus3)ss and (Z0)70(Zr7sCuzs)30 have the highest 
V of 0.372 and 0.373, indicating that these two compositions may have better 
plasticity, which is indeed consistent with the experimental results of three-point 


bending tests. 


3.2 Structural relaxation effect on plasticity 


It is well known that BMGs quenched rapidly from high-temperature melts exhibit 
structural characteristics of relatively low densities and excess free volume. Spaepen 
[25] defines the free volume as the atomic volume in excess of the ideal glassy 
structure described by the densely random packing model. It is generally assumed that 


the amount of free volume affects the plastic flow, and more free volume leads to 
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multiple shear bands and better plasticity [26-28]. For instance, Raghavan ef al. [29] 
observed enhanced plasticity in a Zr-based BMG by increasing the free-volume 
content by ion irradiation, while Murali et al. [30] revealed that a Zr-based BMG 
became brittle with reduced free volume when it was subjected to sub-7, annealing. 
Such a relative change in free volume can be reflected in the enthalpy variation of 
structural relaxation, which has already been proven in previous studies [31,32]. Fig. 
4a shows the DSC curves of (Z0)70(Zr7sCuzs)30 specimens after annealing at 490, 555 
and 620 K below 7, for 1 hour, and the inset is the magnified thermogram describing 
structural relaxation. The data of the as-cast specimen was also shown for reference. 
The as-cast specimen has the largest exothermic enthalpy, suggesting the existence of 
a fairly high amount of free volume. With the increase in annealing temperature, the 
value of exothermic enthalpy gradually decreases, implying that the free-volume 
content drops. In particular, the relaxation peak in the DSC curve of the 620 K- 
annealed specimen becomes invisible. At the same time, the crystallization enthalpies 
remain unchanged, and no crystallization is induced by these annealing treatments. 
The decrease in free volume induced by structural relaxation should also be reflected 
in the relative volume shrinkage [5]. It is found that the density gradually increases 
with the annealing temperature, as shown in Table 2, indicating that structural 
relaxation occurs more sufficiently, which is exactly consistent with the relaxation 
enthalpy change. 

Three-point bending tests were performed to investigate the effect of structural 
relaxation on plasticity. The plastic deflection values decrease obviously with the 
annealing temperature as shown in Fig. 4b, and the specimen becomes totally brittle 
without any visible plastic deformation after annealing at 620 K, in accordance with 
the reduced free volume in the annealed specimens. The ultrasonic measurements 
were further carried out to determine the elastic constants of relaxed 
(Z0)70(Zr7sCuzs)30 specimens, and the values were listed in Table 2. With the increase 
of annealing temperature, v drops from 0.373 for the as-cast state to 0.362 for the 620 
K-annealed specimen, to be also consistent with the deteriorated plasticity. Thus, it 
can be concluded that structural relaxation, which reduces the free-volume content 


and v, results in worse plastic performance in three-point bending test. 


3.3 STZ size determined by nanoindentation creep and statistical distribution analysis 


In order to understand the microscopic mechanisms of various plastic behaviors in 
newly developed Zr-based BMGs and make a thorough inquiry into the correlation 
between STZ size and plasticity, three specimens of Z0, (Z0)70(Zr7sCuzs)30 and 620 K- 
annealed (Z0)70(Zr7sCuzs)30 with different plasticity were selected for the 
nanoindentation tests, and typical indentation load-displacement (P — h) curves were 
shown in Fig. 5. The loading and unloading rates are 5 mN s', and maximum load is 
200 mN with the holding time of 10 s. It is clear that for the as-cast ZO and 
(Z0)70(Zr7sCuzs)30, contact stiffness deduced from the slope of initial unloading data is 
similar, while smaller displacement is reached for Z0 in comparison with 
(Z0);o(Zr;zCuas)ao under the same indentation load, indicating that ZO has higher 
hardness and Young's modulus than (Z0)70(Zr7sCuzs)30. Heat treatment at 620 K does 
not affect the indentation displacement in (Z0)70(Zr7sCuzs)30 obviously, but decreases 
the contact stiffness, suggesting that the contact depth is reduced and the hardness is 
enhanced by annealing. It should be noted that as-cast ZO and 620 K-annealed 
(Z0)70(Zr7sCuzs)30 without bending plasticity had obvious pop-in events, while no 
evident flow serration was observed in the as-cast (Z0)z;o(Zr;sCuas)so with larger 
plasticity, as shown in Fig. 5. Recently, Schuh ef al. [33] found that the decreased 
degree of flow serration in the nanoindentation load-displacement curves acquired on 
a Pd-Ni-P MG was induced by high loading rates or low deformation temperatures, 
consistent with the apparently ductile deformation with more shear bands of smaller 
shear offset at high rates or low temperatures under uniaxial compression tests. Their 
results indicated that the decrease in flow serration was associated with the enhanced 
plasticity, which is further demonstrated by the experimental results of newly 
developed Zr-Cu-Ag-Al-Be BMGs. 

Recently, Pan ef al. [19] characterized the STZs of BMGs by measuring the strain 
rate sensitivity of hardness through rate-jump method based on the Johnson-Sawmer 
CSM model, and the STZ volume can be expressed by 2 = kT/C'mH, where k is the 


Boltzman constant, T is the temperature, m is the strain rate sensitivity, H is the 


1 2RE Gove 1/2 i 
hardness, C = KF em (1 m , R 20.25, € =3, yc = 0.027, Go is the shear 


modulus at 0 K, Tc and ter are the threshold shear resistance at 0 K and temperature 


2/3 
T; = = 0.036, = = | — 0.444 (=) | The strain rate sensitivity m can also be 
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deduced from the creep curve of BMGs here by measuring the dependence of 
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instantaneous hardness on the strain rate [34,35], m = —T = - = ——,, where y 
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is the inelastic shear strain rate, € is the strain rate, o is the critical yield stress. Fig. 6a 
shows the representative displacement-time (h — t) profiles for Z0, (Z0)70(Z1r7sCuzs)30 
and 620 K-annealed (Z0)70(Zr7sCuzs)30 specimens during creep at the maximum load 
of 200 mN, and the fitting curves are described by an empirical law of h(t) = at^ + 
kt, where a, b and k are fitting parameters [36]. It is evident that the empirical law 
can fit the experimental creep data well with high correlation coefficient R? 70.99 for 
all the tests. In comparison with the as-cast and 620 K-annealed (Z0)70(Z1r7sCuzs)30 
specimens, the as-cast Z0 possesses the smaller creep displacement, which agrees 


well with the variation in microhardness listed in Tables 1 and 2. The strain rate of a 


spherical indentation during creep can be approximated as £ = KK [37,38], where 


A is the contact area. The hardness can be obtained by H = P/A = P/7ra”, where the 
contact radius a = J2Rh, — h2 A J2Rh, as he « R, and he = h — CP/S is the 
contact depth at P with C —0.75 for spherical indenter and S the contact stiffness at P. 
The contact stiffness as a function of contact depth during spherical indentation can 
be expressed by S= BE,VA x BEA 2nRh, with ß = 1.129 [39]. Since the 
displacement during creep is much smaller than that during loading prior to creep, the 
contact stiffness during creep can be roughly considered as a constant. Fig. 6b shows 
In H as a function of In & during creep for Z0, (Z0)70(Zr7sCuzs)39 and 620 K-annealed 
(Z0)70(Zr7sCuzs)30 specimens, respectively. From the slope of the curves in Fig. 6b, 
the value of strain rate sensitivity m can be obtained. In crystalline materials, only the 
creep data in the steady-state stage was adopted, and people utilized the dependence 
of hardness on steady-state creep rate to calculate the strain rate sensitivity and 
described the creep deformation kinetics [40]. Following this empirical rule, only 2 
calculated from the steady-state creep stage of BMGs was discussed. In Fig. 6a, the 
creep deformation becomes steady with almost constant creep rate in the time region 
of 6-10 s for each specimen. Therefore, the typical m calculated from the linear fitting 
of the steady-state creep stage in Fig. 6b is listed in Table 3, and the corresponding 2 
can be obtained according to N = KT /C'mH, where the hardness at the end of creep 
process is used for H due to its slight change during steady-stage creep. It is clear that 
the value of m in (Z0);o(Zr;zsCuos)so decreases as compared with that in ZO due to 


minor alloying, and annealing treatment at 623 K further reduces m of 
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(Z0)70(Zr7sCuzs)30. The calculated 2 rises from 5.51 nm? for Z0 to 7.25 nm? for 
(Z0)70(Zr7sCuzs)30, and further increases to 9.00 nm? in 620 K-annealed 
(Z0)70(Zr7sCuzs)30. Once the value of is determined, the number of atoms involved 
in a STZ event can be subsequently assessed by the dense-packing hard-sphere model 
for MGs. The average atomic radius is statistically estimated as r = ($? A; 7) /? 
[19], in which A; and r; are the atomic fraction and the atomic radius of each element, 
respectively. The variation in the representative number of atoms in STZs with minor 
alloying and annealing treatment is also listed in Table 3. 

Since more than 400 indentations on the surface were conducted for each 
composition, the distribution of the creep displacement and hardness indeed exists. 
Consequently, the value of 2 deduced from each nanoindentation test varies from 
each other, and the distribution of) for each specimen can be obtained as shown in 
Fig. 6c. The peak position gradually shifts from 5.50 nm? for Z0 to 7.50 nm? and 8.70 
nm? for the as-cast and 620 K-annealed (Z0)70(Zr7sCurs)30, respectively, which is 
consistent with the data from the typical creep curves listed in Table 3. It is obvious 
that from the as-cast Z0 to (Z0)70(Zr7sCuzs)30, the average value of 2 increases and is 
positively correlated with plasticity and v, according with the proposition of Pan ef al. 
[19], 1.e., the ductile glass has the large 2. However, after annealing treatment at 620 
K, (Z0)70(Zr7sCuzs)30 becomes brittle with a lower v, but the value of 2 is further 
enlarged. 

Recently, Schuh ef al. [41] investigated the statistical distribution of rate-dependent 
pop-in events in the nanoindentation curves of SiC material, and associated the first 
pop-in with the onset of plastic deformation and the defect nucleation. In the 
framework of stress-assisted, thermally activated defect nucleation, they proposed a 
nucleation theory-based analysis for incipient plasticity during nanoindentation, and 
deduced the activation volume of defect nucleation events. Since the STZ event of 
MGs can also be considered as a stress-biased thermally activated process, Choi ef al. 
[38] studied the cumulative distribution of the maximum shear stress at first pop-in 
Tmax by adopting the Schuh's method, and the value of can be obtained. In 
spherical indentation, the maximum shear stress at the elastic-to-plastic transition, 


Tmax © 0.45P,,, appears at a distance of half the contact radius below the center of 
1 
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lrir = In the work of 
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the indenter [42], where the mean contact pressure Pn = ( 


12 


Choi et al. [38], Tmax at the first pop-in was considered as the critical shear strength 
(7y) for the onset of plasticity in the indented materials. However, Packard and Schuh 
[43] recently pointed out that (1) high local stress at a local point in MGs cannot 
initiate a shear band, and shear banding occurs only when the yield strength is 
exceeded everywhere along a potential shear band path; (ii) previous analysis to 
extract shear yield stress from Hertzian contact theory overestimate the shear yield 
strength, and ty = 0.07F,,. The scientific debate still remains over the exact relation 
between T, and Pn, but ty = H/ 3v3 [19] in MGs is generally accepted. In order to 
confirm the variation in 2 with minor alloying and annealing treatment mentioned 
above, we calculated Q by analyzing the cumulative distribution of t, deduced 
from H instead of Pn. For each composition, more than 400 nanoindentation tests 


were performed in order to obtain the statistical distribution of ty. In the spherical 


: : . >  dP d(H Pd 
nanoindentation, the loading rate P = — = gon = 2H na“ <-2-Z , and 
dt dt a dt 
1 dVA ` 1da P 
¿= —=— = Therefore, the equivalent strain rate can be expressed by £ = — 


VA dt adt 2P' 
and is equal to 1.25x10? sl. The loading-unloading curves were used to calculate 
hardness values with the histogram of the hardness distribution shown in Fig. 7a. The 
change of average hardness with composition and annealing treatment is consistent 
with the results of microhardness tested by the Vickers indenter. The corresponding 
cumulative distribution of ty ~ H/ 3V3 for Z0, (Z0)70(Zr7sCuzs)30 and 620 K- 
annealed (Z0)70(Z1r7sCuzs)30 was illustrated in Fig. 7b, respectively. Schuh et al. [41] 
suggested that the thermally-assisted and stress-biased yielding always exhibits a 
spread in the yield strength due to the fact that thermal noise sometimes favors 
yielding and sometimes works against it. On this basis, the cumulative probability of 


elastic-to-plastic transition can be described as a function of shear stress f = 1 — 


k -AF* V" h [i 
bê [- TM exp ( — ) exp EA where N, is the attempt frequency, 7 is the stress 


Zë 


rate, V* and AF" are the activation volume and energy of the STZ event. It can be 
rewritten as In[In(1 — f) ! I=Ệ + In (— 5) + = and the value of V* can be 


obtained from the slope of In[In(1 — f)~*] versus ty, as shown in Fig. 7c. Only the 
data of f ranging from 20% to 80% was adopted to ensure the availability of each 
linear fitting. Based on the CSM model, the STZ volume can be expressed by 
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Tc 


= Rare eV [19]. The values of V* and 2 of various specimens were 
listed in Table 3. It is clear that the variation tendency of 2 after minor alloying and 
annealing treatment is well consistent with that observed in the steady-state creep 
experiments. Meanwhile, as compared with .2 from the steady-state creep with lower 
strain rate (5x10'~1x103 s^), the value of Q is slightly decreased due to the higher 
strain rate (1.25x10? s) adopted in the cumulative distribution method. Such a 
dependence of 2 on the strain rate was already observed by Choi er al. [38], who 
investigated statistical distribution of the first pop-in stress in a Zr-based BMG at 
different loading rates and found that increasing the strain rate can obviously reduce 
N, although some inconsistent results on the influence of strain rate on Q were 
reported [44]. 

As for the effect of structural relaxation on the plasticity, there was almost no 
debate upon the deteriorated plasticity of the relaxed MG [28], while some 
controversial results on the influence of structural relaxation on Q existed. For 
example, Pan et al. [45] and Jiang et al. [16] found that 2 decreased with annealing, 
and the annealed MG specimen became brittle. During structural relaxation, the free- 
volume content was decreased, indicating the possible sites available for STZ 
nucleation reduced. Meanwhile, since 2 decreased due to structural relaxation, more 
STZs were required to be activated to reach the percolation limit for shear banding, 
suggesting that the plasticity of relaxed specimen would deteriorate. Nevertheless, 
Choi et al. [46] revealed that annealing treatment augmented 2 since a more uniform 
structure might enable larger number of atoms to be involved in the cooperative 
shearing for a STZ event. Although 2 increased with annealing, the reduced 
nucleation sites for STZs due to the decreased free volume would eventually result in 
the brittle fracture during deformation. Based on the results mentioned above, 2 may 
be not the only factor controlling the plasticity of MGs, and free volume providing the 
possible sites for STZ nucleation can also play an important role in affecting the 


plastic deformation behavior. 


3.4 Packing efficiency from SR-XRD and free volume from PALS 


To further understand the enhanced plasticity through tuning composition and the 


reduced plasticity through annealing treatment from a structural aspect, SR-XRD 
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experiments were performed. Fig. 8 shows the structural factors S(q) for ZO, 
(Z20)80(2r75Cu25)20, (20)70(Zr75Cu25);o and 620 K-annealed (Z0)70(2r75Cu25);0. With 
the increase in the fraction of Zr;sCuas, the first peak of S(q) shifts gradually towards 
a lower q value, indicating an increase in average atomic distance due to the higher 
content of Zr with the largest atomic radius among the constituent components, but 
the peak profile does not change significantly. And the difference in S(q) between the 
as-cast and 620 K-annealed (Z0);o(Zr;zCuas)ao specimens is almost negligible. 
Considering the peak in the nearest atomic pair, reflecting the short-range order, the 
magnified curves of the first peak of G(r) are shown in Fig. 9a. Evidently, in 
comparison with Z0 and (Z0)so(ZrzsCuos)oo, the first peak of (Z0)70(Zr7sCuzs)30 shifts 
to a higher r value. Using the dense-packing hard-sphere model, the interatomic 
distances and concentration-weighted partial radial distribution values for ZO, 
(Z0)so(Zr;szCuas)ao and (Z0)79(Zr7sCuzs)30 are indicated in Fig. 9b-d, and also listed in 
Table 4. It is clear that Zr-Zr, Zr-Cu and Zr-(Ag, Al) atomic pairs dominate the first 
peak. The relative proportion of the Zr-Zr atomic pair with the largest bond length in 
(Z0)70(Zr7sCuzs)30 specimen is obviously higher than those in other two compositions 
owing to its higher Zr content, leading to the shift of the first peak of G(r) to a higher 
r value. The coordination number (CN), a key parameter to evaluate the atomic 
packing efficiency, can be deduced by integrating G(r) curve over the range of the 
first peak. The calculated CN was 14.9, 13.3, 13.2 for Z0, (Z0)so(Zr;zCuas)ao and 
(Z0)70(Zr7sCuzs)30, respectively. 

Combining the average atomic distance with the value of CN, it can be intuitively 
understood: with the increase of Zr content, packing efficiency drops. Such a drop in 
atomic packing density is induced by composition tuning, mainly due to the increase 
in the amount of Zr with the largest atomic radius. Actually, it does not indicate for 
sure that the free-volume content increases with the fraction of Zr;zCuas. At the same 
time, the difference in atomic packing between the as-cast and 620 K-annealed 
(Z0)70(Zr7sCuzs)30 specimens is not obvious from the SR-XRD observation, but 
annealing-induced embrittlement indeed occurs. Generally, structural relaxation 
reduces free volume in MG, suggesting that S(q) and G(r) curves are not enough 
sensitive to detect the variation of free volume caused by annealing. Therefore, PALS 
measurements were carried out to gain some information about free volume. 


Experimental positron lifetime spectra obtained in Z0, (Z0)so(Zr75Cu25)20, 
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(Z0);o(Zr;zCuas)ao and 620 K-annealed (Z0);o(Zr;zCu;s)ao are shown in Fig. 10. 
Analysis of annihilation spectra reveals that two lifetime components fitting to the 
data is the best possible fit for each specimen, corresponding to two positron trap 
sizes. Since Flores ef al. [47] recently identified three characteristic lifetimes for a Zr- 
Cu-Ni-Al-Nb BMG, we also tried to fit the data by using three lifetime components, 
but such fit cannot satisfactorily describe the positron annihilation spectra with the 
variance of the fit greater than 2. With the increase in the fraction of Zr7sCuzs in the 
as-cast state, the average lifetime for the shorter lifetime component (7¡) gradually 
decreases from about 187 ps for Z0 to 166 ps for (Z0)s0(Zr7sCuzs)20 and 164 ps for 
(Z0)70(Zr7sCuzs)30, as listed in Table 5. Meanwhile, the normalized intensity of T,, 
indicative of the relative concentration of positron trap, increases from 91.6% for ZO 
to 94.9% for (Z0)go(Zr7sCu2s)20 and 95.7% for (Z0);o(Zr;sCu;s)so. These results 
suggest that the size of smaller open volume decreases, but its concentration increases. 
The shorter lifetime 7¡ is similar to the lifetime in bulk states of crystalline Zr, and 
can be attributed to annihilation at tetrahedral and octahedral holes in the densely 
packed MGs [48]. In the Zr-Cu-Ni-Al-Nb BMG, Flores et al. [47] found that the open 
volume with the lifetime (~0.3 ns) comparable to that of a monovacancy in crystalline 
Zr existed, and these regions were referred to as flow defects. In the BMG specimens 
investigated here, no such lifetime component was detected, but it did not mean that 
flow defects were absent, otherwise plastic deformation should not be observed. In the 
densely-packed glass structure, tetrahedra, possessing highest packing efficiency, 
cannot be combined to fill space, and nontetrahedral Bernal-type polyhedra are 
geometrically necessary. These tetrahedral and nontetrahedral holes with size much 
less than that of a vacancy are not flow defects according to the classification of 
Flores et al. [47], but can be considered as flow defect nuclei, which may lead to the 
formation of new flow defects with additional distortion under stress and act as 
initiation sites for STZs. The longer lifetime component (725) is of sufficient 
magnitude to be associated with the annihilation of ortho-Positronium (o-Ps) in voids. 
With the increase in the fraction of Zr7sCuzs, T, gradually rises from about 2.009 ns 
for Z0 to 2.481 ns for (Z0)s0(Zr7sCuzs)20 and 2.741 ns for (Z0)70(Zr7sCuzs)30, as listed 
in Table 5. The normalized intensity of r; drops from 8.4% for Z0 to 5.1% for 
(Z0)so(Zr;sCuas)a2o and 4.3% for (Z0)70(Zr7sCuzs)30. The average lifetime for the whole 
sample gradually reduces from 0.339 ns for Z0 to 0.283 ns for (Z0)g0(Zr7sCuzs)20 and 
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0.274 ns for (Z0);o(Zr75Cuos)so. Since the lifetime of o-Ps is directly related with the 
size of the voids in which it forms [49], the radius of the void is calculated to be 2.86, 
3.26 and 3.46 A for Z0, (Z0)g0(Zr7sCuzs)20 and (Z0)70(Zr7sCuzs)30, respectively, much 
larger than that (~1.6 A) of Zr atom. Such large voids are possibly not involved in the 
STZs nucleation [17] since its radius is almost one third of the STZ radius obtained 
experimentally. The existence of voids may mainly trigger the catastrophic rapid 
propagation of shear bands and consequently deteriorate the plasticity. 

In (Z0)70(Z175Cuzs)30, shorter lifetime component (7¡) decreases from 164 ps for the 
as-cast state to 156 ps for 620 K-annealed state, while no obvious changes in the 
longer lifetime component (T2) and the normalized intensities of 7¡ and v; are 
detectable, as listed in Table 5. It indicates that average lifetime is decreased by 
structural relaxation and positively correlated with the plasticity. Note that from Z0 to 
(Z0)70(Zr7sCuzs)30, the bending plasticity rises, while the average lifetime drops due to 
the reduction in average free volume, suggesting that the average lifetime (average 
free volume) is not always correlated positively with the plasticity. Meanwhile, with 
the increase in the fraction of Zr;sCuas, the plasticity is enhanced and 2 becomes 
larger, while structural relaxation results in the embrittlement accompanied with the 
further enlarged 2. The question is why these two conflicting trends of plasticity with 
average lifetime or exist, and rationalizing the evolution of plasticity with the 
composition and structural relaxation from the combined viewpoint of open volume 
size distribution and 2 is of significance for us to understand the microscopic 
mechanisms of plastic deformation in MGs. It is proposed that the smaller open 
volume is responsible for the nucleation of STZ event, while the larger open volume 
(voids) possibly induces the catastrophic propagation of shear band (cooperative 
movement of STZs). From Z0 to (Z0)70(Zr7sCuzs)30, increases by 1.316~1.340 
times, as listed in Table 3, and the concentration of smaller open volume increases by 
1.04 times. Supposing the distribution of open volume in MGs is uniform, one STZ 
event in (Z0)70(Zr75Cu25)30 Consumes more number of smaller open volume as 
compared with Z0. Because the increase in the concentration of smaller open volume 
cannot compensate it, the opportunity for STZs to reach the percolation limit for shear 
banding is reduced and subsequently the plasticity becomes worse, which is in 
contradiction to the experimental results. Thanks to the larger open volume regions, 


the situation may be changed conversely. The reduced concentration of larger open 
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volume from ZO to (Z0)70(Zr7sCuzs)30 suggests that the chance for shear bands to meet 
the voids (1.e., larger open volume) during their propagation lessens. Meanwhile, the 
contribution of larger open volume to the average positron lifetime is about 169 ps for 
Z0 and 118 ps for (Z0)70(Zr7sCuzs)30, respectively. Therefore, the decrease in the large 
voids in (Z0)70(Zr7sCuzs)30 inhibit the catastrophic runaway of shear bands and 
enhance the plasticity accordingly, although the size of voids in (Z0)70(Z1r7sCuzs)30 
increases as compared with that in Z0. The influence of larger open volume may be 
predominant during the competition with the influence of smaller open volume, 
eventually resulting in the enhanced plasticity in (Z0)70(Zr;5Cu»5)so as compared with 
ZO. Structural relaxation at 620 K for 1 hour in (Z0)zo(Zr;zCuas)so drops only the size 
of smaller open volume without any change in the concentration, but increases 2 by 
1.241~1.417 times, as listed in Table 3. As mentioned above, the increased AU 
consumes more number of smaller open volume, while the concentration of smaller 
open volume does not vary with the annealing treatment. Therefore, there is less 
chance for STZs to reach the percolation limit for shear banding in the 620 K- 
annealed specimen, which becomes brittle with the reduced average free volume, 
although it possesses larger 2 than the as-cast specimen. Based on all the results 
shown here, the open volume size distribution may be more crucial than the average 
lifetime for the plastic deformability of BMGs, and considering both 2 and open 
volume size distribution simultaneously can explain the evolution of plasticity with 


the composition and structural relaxation. 
4. Conclusions 


In summary, a comprehensive investigation into the effect of compositional tuning 
and structural relaxation annealing on the microstructure and mechanical properties in 
Zr-Cu-Ag-Al-Be BMG was performed. The results allowed the following conclusions 
to be drawn: 

(1) The bending plasticity of the as-cast BMG is gradually enhanced with increasing v 
and decreasing GFA. Meanwhile, structural relaxation worsens the plasticity with the 
reduced free volume and v. All the results indicate that the plasticity correlates closely 
with v, which is a good apparent indicator of plasticity. 

(2) Through the creep method and statistical analysis method, (2 was determined for 


various BMGs. With the increase in v and plasticity from Z0 to (Z0)70(Z1r7sCuzs)30, the 
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value of 2 determined from both methods is found to be increased. However, after 
annealing treatment at 620K, (Z0)70(Zr7sCuzs)30 becomes brittle with a lower v, but 
the value of 2 is further enlarged. All the results suggest that only 2 cannot 
satisfactorily describe the plastic deformability of BMGs. 

(3) By the PALS measurements in Z0, (Z0)so(Zr;5Cu»s5);o and (Z0)70(Z17sCu2s)30, two 
lifetime components, corresponding to two positron trap sizes, provide the best 
possible fit for each specimen. It is revealed that with the increase in the fraction of 
Zr75Cu»s, the size of smaller open volume (tetrahedral and octahedral holes) decreases 
with increased concentration, and the size of larger open volume (voids) increases 
with decreased concentration. Both the contributions of smaller and larger open 
volume to the average positron lifetime drop with the increase in the fraction of 
Zr75Cuos. Structural relaxation at 620 K in (Z0)0(Zr75Cu25)30 only reduces the size of 
smaller open volume without changing the concentration, leading to the decreased 
average positron lifetime. All the results indicate that average positron lifetime does 
not always correlate positively with the plasticity. 

(4) From ZO to (Z0)70(Zr7sCuzs)30, more number of smaller open volume is occupied 
for the STZ activation due to the increased 2, and the rise in the concentration cannot 
compensate it, leading to the reduced opportunity for STZs to reach the percolation 
limit for shear banding and consequently deteriorated plasticity. On the other hand, 
the concentration of large voids and its contribution to the average positron lifetime 
also decrease with the increase in the fraction of Zr;sCuos, and the catastrophic 
runaway of shear band is arrested, enhancing the plasticity accordingly. The influence 
of the larger voids on the propagation of cooperative STZs may be predominant over 
the influence of smaller open volume on the nucleation of STZs, eventually resulting 
in the enhanced plasticity with the increase in the fraction of Zr7sCuzs. Since structural 
relaxation at 620 K in (Z0)+o(Zr;zCuas)so only reduces the size of smaller open volume, 
more number of smaller open volume is required for STZs activation due to the 
increased (2, and there is less chance for STZs to reach the percolation limit for shear 
banding. Thus, the (Z0)70(Zr7sCuzs)30 specimen annealed 620 K becomes brittle with 
reduced open volume, although it possesses larger than the as-cast specimen. AII 
the results indicate that both Q and open volume size distribution play crucial roles in 


affecting the plasticity of BMGs. 
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Table captions 
Table 1 
Thermal, mechanical properties and density parameters of the as-cast specimens with 


Z0 representing Zr46Cu3014A gs 3c AlgBe; 5. 


Table 2 
Mechanical properties and density changes of (Z0)70(Zr7sCuzs)30 after annealing for 1 


hour at different temperatures. 


Table 3 
STZ size calculated based on the creep data and the statistical analysis of the yield 


data, respectively. 


Table 4 
Interatomic distances (D) estimated from the hard sphere model and concentration- 


weighted partial radial distribution values (W;;) for three as-cast specimens. 


Table 5 


Positron annihilation lifetime spectra (PALS) for the as-cast and annealed specimens. 
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Figure captions 
Fig. 1.  Pseudo-ternary composition map with Z0 representing 


Zr46Cu3z0.14Ags36AlgBe7 5. 


Fig. 2. XRD patterns (a) and DSC curves (b) of the as-cast plate specimens with the 
thickness of 1.5 mm. 


Fig. 3. Flexural stress vs. deflection curves (a) of the as-cast plate speicmens with 
various compositions; SEM images of tensile surface after fracture for (b) Z0, (c) 


(Z0)ssZr¡s, (d) (Z0)es(Zre;7Cus3)ss and (e) (Z0)7o(ZrzsCuos)so, respectively. 


Fig. 4. DSC curves (a) and flexural stress vs. deflection curves (b) of 
(Z0);o(Zr;szCuas)ao plates after annealing at different temperatures for 1 hour. The data 


of the as-cast specimen is shown for reference. 


Fig. 5. Room-temperature nanoindentation load vs. displacement curves for the as- 


cast and the annealed specimens. The curves are relatively shifted for clarity. 


Fig. 6. Typical experimental and fitting curves (a) during creep at 200 mN load; The 
relationship (b) between hardness and strain rate during creep; The distribution of 
STZ (c) obtained by nanoindentation creep tests. The circles marked in (a) and (b) 


indicate the steady-state creep stage. 


Fig. 7. The distribution of hardness (a), cumulative probability distribution of shear 


yield stress (b), and In[In(1-f) !] vs. shear yield stress (c) for various specimens. 

Fig. 8. Structural factors S(q) for various specimens. 

Fig. 9. Magnified curves (a) of the first peak of G(r) for various as-cast specimens; 
(b)-(d) the first peak of G(r) for three specimens with the bar indicating the weight 


values of atomic pairs. 


Fig. 10. Experimental PALS curves obtained in various specimens. 
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